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SUMMARY 


The incompressible three-dimensional, turbulent boundary layer (3DTBL) experi- 
ments of van den Berg and Elsenaar (ref. 1), Dechow (ref. 2) and MUller (ref. 3) 
were simulated numerically by integrating the boundary-layer equations together with 
an algebraic eddy-viscosity turbulence model. For the flow treated by van den Berg 
and Elsenaar, the downstream portion, where the crossflow was large, could not be 
predicted with the present computational method; we feel that this flow was signifi- 
cantly influenced by elliptic flow-field effects. Though Dechow's experiment also 
indicated departures from the boundary-layer concept, our calculations agreed reason- 
ably well with the mean-flow development up to separation. In MUller' s experiment 
the normal pressure gradients were found to be negligible in regions with large skew- 
ing and enabled us to test turbulence models using the boundary-layer equations. 

The simulation of this flow compared favorably with the experimental data throughout 
the flow field and suggested the applicability of algebraic eddy-viscosity models for 
3DTBLs. 


INTRODUCTION 


Predictions of 3DTBLs are needed for many engineering purposes, for example, in 
designing swept, low aspect ratio wings and improving their performance. However, 
the development of computational methods, as well as our understanding of the turbu- 
lent momentum transfer in three-dimensional flows, has lagged substantially behind 
the state of art reached for two-dimensional flows. Since only a small number of 
experiments was available, which documented both the mean flow field and the Reynolds 
stresses of 3DTBLs, turbulence modeling for such flows was not always confirmed 
experimentally. For most numerical computations of 3DTBLs, turbulence models devel- 
oped for two-dimensional flows have been applied together with an ad hoc closure 
assumption for the crossflow momentum equation. The predictions of low speed flows 
with moderate crossflow generally compared well with experimental data, but flow 
fields with large crossflow often could not be simulated with reasonable accuracy 
(ref. 4), especially the flow field experiment performed by van den Berg and 
Elsenaar (ref. 1). Comparison computations for this experiment failed to predict the 
three-dimensional part of the flow ahead of separation. Several attempts have been 
made to improve the results of this experiment. Krause (ref. 5) and Kordulla, for 
example, investigated the applicability of several algebraic turbulence models, but 
found the results of the computations to be strongly coupled with the numerical 
accuracy and also with the accuracy with which the outer-edge boundary conditions 
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were prescribed. Since isotropic turbulence models assume equal directions of the 
resultant, turbulent shear stress vector and the rate-of-strain vector throughout the 
flow field, Rotta (ref. 6) used the measurements of Elsenaar and Boelsma (ref. 7), to 
develop an anisotropic closure assumption for the pressure-strain correlation and 
obtained slight improvement of the predictions. From the same set of measurements 
Cousteix et al. (ref. 8) developed a heuristic transport equation for the direction 
of the shear stress vector; their computations also indicated some improvements of 
the results only. Rubesin's computations as discussed in reference 9, used a closure 
model for the transport equations of the Reynolds stresses and did not require an \ 
a priori assumption for the degree of anisotropy, but again for the three-dimensional 
flow region only small improvements were achieved compared to the other calcula- 
tions. All predictions mentioned so far did not only underestimate the development 
of the secondary flow, but also failed to predict separation. The latter could only 
be enforced under modified conditions like those of Elsenaar et al. (ref. 10), who 
increased the prescribed pressure distribution and also reduced the influence of the 
turbulence model by adjusting several coefficients to the experimental data. 

Based on the results of 3DTBL simulations like those mentioned above, Marvin 
(ref. 9) suggested that eddy-viscosity turbulence models may not be adequate for com- 
puting flows with severe crossflow but need considerable improvement. The present 
author feels that further comparisons of measured and predicted mean-flow and turbu- 
lence fields are needed before discarding simple, gradient-diffusion turbulence 
models. In fact, during the recent years several documentations of 3DTBL develop- 
ments including the measurements of all components of the Reynolds stress tensor were 
completed, for example by Dechow (ref. 2) and Muller (ref. 3). To the author’s 
knowledge these experiments have not yet been used for comparing theory and experi- 
ment . 


In the present study the three flow fields were simulated numerically by solv- 
ing the boundary-layer equations. We concentrated our efforts on investigating 
whether an algebraic turbulence model as developed for two-dimensional flows allows 
for reasonable predictions of 3DTBLs. In order to sort out the applicability of a 
particular turbulence model, we also had to study the sensitivity of the predictions 
to modifications of the closure assumptions and to the accuracy of the pressure dis- 
tribution prescribed. 


NUMERICAL METHOD AND TURBULENCE MODELS 


The equations for incompressible 3DTBLs were integrated in a Cartesian frame of 
reference: 

3U/3x + 3V/3y + 3W/3z = 0 ; 3p/3y = 0 


U3U/3x 4- V3U/3y + W3U/3z = -(l/p)3p/3x 

+ 3 (v3U/3y - uv)/3y 
U3W/3x + V3W/3y + W3W/3z = -(l/p)3p/3z 

+ 3(v3W/3y - vw)/3y 


( 1 ) 
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Quasi-two-dimensional flows were computed setting 3/3z = 0; while for two-dimensional 
flows W was set to zero. The boundary conditions included no-slip conditions at the 
wall and prescribed outer-edge velocities, which were determined by smoothing the 
measured pressure distribution and solving the Euler equations. The baseline turbu- 
lence model used for closing equation (1) was formed by Michel et al.’s (ref. 11) 
mixing-length formula together with the assumption of isotropic eddy viscosities. 

-uv - v 9U/3y ; -vw = v 9W/9y ; v = 

X z c 

v = 4 2 F 2 {(aU/3y) 2 + (aw/ay ) 2 }' 5 ; F = 


4/6 = 4 e /6*tanh{0.41/(4 e /6)*(y/6)} ; u t 

We applied Cebeci's (ref. 12) pressure gradient correction to the van Driest 
damper F with s being the coordinate along the outer-edge streamline: 

N = (1 - 11.8 p + )^ ; p + = -{1/ (pu^) }ap/8s . (3) 


v = v 
X z 


1 - exp{-yu N/(26v)} 

T 


\J/p s W 6 = °- 085 J 


( 2 ) 


Additionally, we used a closure assumption with anisotropic eddy-viscosities. As 
suggested by Rotta (ref. 6) the Reynolds stresses were defined in a local coordinate 
system, the x-axis of which was aligned with the yaw direction of the resultant flow 
vector. Then the stresses were transformed into the frame of reference used for 
equation (1) . 


-u v = v 9U /3y 
mm x m 


-v w 
m m 


Tv z aw ra /a y 


= v, 


v = Tv 
z t 


v = 4 2 F 2 {(3U /3y) 2 + T(8W /8y) 2 } 
t m m 


H. 


(4) 


The ratio T of the eddy viscosities defined the degree of anisotropy. 


The initial conditions for the spanwise integration were obtained either from 
solving the plane-of-symmetry equations (ref. 12) or by integrating equation (1) 
with Hall’s difference scheme. The CFL-condition was applied in the wall-streamline 
direction. Then the integration was carried out with Crank-Nicholson differencing. 
For all computations we used step sizes Ax = Az = 8mm; the normal grid spacing Ay 
was smaller than y+ = 1 at the wall and then was stretched towards the outer edge 
yielding 70 to 150 node points across the boundary layer. As found from test calcu- 
lations, grid effects were negligible for all results obtained. 


The integration of the governing equations was- based on the algorithm of Krause 
et al. (ref. 13). Using the baseline turbulence model the computational method was 
tested against case A4 of reference 4 (turbulent boundary layer approaching a cylin- 
der mounted perpendicular to the surface). As shown in figure 1, the results 
obtained for the turning angle 3 W between outer-edge and wall streamline directions, 
and for the shape factor H evaluated with the streamwise velocity profile, agreed 
well with the predictions of other differential methods. 
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RESULTS OF NUMERICAL SIMULATIONS 


Experiment of van den Berg et al. (Refs. 1, 7) 

This investigation of a 3DTBL on a swept flat plate was performed under 
"infinite swept-wing conditions", that is vanishing gradients of all flow variables 
in spanwise direction, figure 2. By means of a diffuser-type channel an adverse 
pressure gradient normal to the leading edge was generated, such that three- 
dimensional separation occurred in the downstream region. 

We obtained the initial conditions for starting the numerical simulations by 
matching a two-dimensional, zero-pressure gradient, turbulent boundary-layer calcula- 
tion with Reynolds number Re =2.42 x 10 6 (U^ defined at station 1; L = lm) to the 
measured momentum-thickness Reynolds number and skin friction coefficient. Then a 
small crossflow yielding B w 1.5 was prescribed according to reference 1. 

The computation of the downstream flow development was carried out under quasi- 
two-dimensional conditions with 9/8z = 0. Using the pressure distribution measured 
at the wall together with the baseline turbulence model we obtained the result indi- 
cated by (run Al) in figures 3 and 4. Up to station 5 the measured and calculated 
mean velocity profiles as well as the skin friction coefficients C£ = 2x w /(pU|) 
and the wall turning angles (a+3 w ) agreed well. The angle a defined the direction 
of the outer-edge flow with respect to the x-axis. Further downstream the compari- 
son indicated deviations which increased rapidly; at station 7 the computed wall 
turning angle lagged 15° behind the measured value of 40° , while the mean velocity 
component_U was overpredicted as much as 100% in the near-wall region. The spanwise 
velocity W was obtained a few percent too large. Beyond station 7 the computed 
values of Cf and $ w were approximately constant and did not lead to separation. 

These results were identical to those reported earlier, for example by Krause (ref. 

5) . Carrying out a fully three-dimensional calculation by using the corresponding 
wall pressure distribution of reference (1), we did not observe any improvements in 
the predicted results. When using a standard turbulence model the calculation was 
insensitive to the departures from the ideal infinite swept wing case. 


In the three-dimensional flow at station 7 the Reynolds shear stresses u m v m 
and v m w m acting in the local yaw direction and normal to it, respectively, indicated 
peak values at 12 mm < y < 15 mm, which were 25 to 30% higher than those of the 
experimental data. Since in the outer layer the calculated mean velocity gradients 
were evaluated smaller than the measured ones, these results seemed to indicate an 
overestimation of the turbulent viscosities by the closure assumption used. That is 
why we reduced the magnitude of viscous forces in subsequent calculations. First we 
considered Rotta's definition of anisotropic eddy viscosities, equations (4). In 
reference 7 the ratio of the local crosswise to streamwise viscosity was measured as 
0.5 < T < 0.6 upstream of station 5 and as T > 0.7 downstream of this station. 
According to Mliller’s (ref. 14) error analysis for turbulence measurements in 3DTBLs, 
accurate results for this ratio are very difficult to obtain, especially in weak 
three-dimensional flows. That is why for the present calculations the ratio T was 
set equal to one initially and then was gradually reduced beyond station 4 or 5, 
respectively, such that the prescribed experimental value was obtained at the next 
downstream measuring station. In our calculation with T = 0.7, introduced at sta- 
tion 5, we found negligible differences compared. to run Al. 
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Cousteix et al. (ref. 8) applied T < 1 throughout the flow field; their com- 
putations yielded, for example, overestimated turning angles in the flow region with 
moderate crossflow and then seemed to improve the results slightly in the three- 
dimensional part. 

Second, we investigated the sensitivity of the results to the mixing-length 
distribution prescribed. While in reference 7 the measured inner-layer lengths 
scales were close to i = 0.41y, the normalized outer-layer data approached 
Z /6 = 0.07 at station 5 and 0.05 at station 7; thereby they deviated significantly 
from the distribution assumed by the baseline model. Plotting the mixing length 
distributions £(y), however, indicated that these profiles were approximately 
constant downstream of station 4 or 5 and did not scale with the boundary-layer 
thickness 6. "Freezing" the length scale distribution at one of these stations 
effectively was_identical to the empirical relationship of Elsenaar et al., reference 
10, £ e /6 S (1-V e ) where V e is the normal velocity at the outer edge. When using 
equation (2) up to station 5 and then "freezing" Z(y ) , the results were close to 
those of run A2, figure 3. For run A2 T = 0.7 was used in addition to the frozen 
mixing length. Downstream of station 6 the skin friction coefficient decreased 
slightly with increasing x-coordinate, while between stations 7 and 8 B w was 5 
higher than in run Al . Again, separation was not predicted. 

In further computations we investigated the sensitivity of the predictions to 
the prescribed static pressure distribution. In reference 1 the pressure at the 
outer edge of the boundary-layer was derived from the mean-flow measurement and was 
found to be larger than that at the wall downstream of station 6. For our test cal- 
culations we assumed a pressure rise between stations 6 and 7 leading to a value that 
was 2% higher than the one measured at the wall. The downstream outer-edge pressure 
data were approximated by a polynomial curve fit. The approximation of the outer- 
edge data as indicated by the broken line in figure 2 will be used throughout the 
calculations discussed hereafter. Using the baseline turbulence model the results 
indicated an increase in $ w by 4° , while Cf stayed about the same. However, when we 
applied the larger pressure distribution together with the reduced eddy viscosities 
of run A2, we obtained a monotonic decrease of Cf (run A3 in figure 3) and separa- 
tion occurred between stations 9 and 10. At station 7 the spanwise velocity W agreed 
well with the experimental data, figure 4, while the velocity component U was still 
overpredicted. The turbulent shear stress profiles decreased considerably compared 
to the previous runs; except for a narrow region close to the wall they were smaller 
than in the experiment. Introducing T = 0.7 at station 4 had little effect com- 
pared to run A3 where T = 0.7 was introduced at station 5. 

For run A4 we used isotropic eddy viscosities and kept the mixing length dis- 
tribution constant downstream of station 4; then we also introduced T = 0.7 at 
the same station, run A5. Compared to previous computations Cf decreased further, 
and B w increased substantially. For run A5 B w even indicated an overshoot compared 
to the measurements, while both mean velocity components were close to the experi- 
mental data. For this case the momentum thickness, defined with the streamwise 
velocity profile in the direction of the resultant outer-edge velocity, was evalu- 
ated at station 7 as 5.95 mm compared to the experimental value of 6.35 mm (5.06 mm 
in run Al) , while the shape factor H = 1.65 was identical to the experimental 
value (1.45 in run Al) . At the s ame s tation the peaks in the profiles of both 
Reynolds shear stresses u m v m and v m w m were about 30 percent below the measured ones. 
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Dechow's Experiment (Ref. 2) 


Dechow investigated a 3DTBL approaching a tear-drop shaped body extending 
between lower wall (test surface) and upper wall of a channel with a height of 
300 mm, figure 5. Most measurements of mean velocities and Reynolds stresses were 
taken at stations oriented along an external streamline. Since the present computa- 
tional method was limited to attached flows, we could not simulate any measurements 
downstream of station 6. The two-dimensional initial conditions at x = 0 were 
obtained from smoothing the measured velocity profile U(y) of station 1. We matched 
the momentum thickness Reynolds number and the skin friction coefficient to the 
experimental data. The free-stream Reynolds number was prescribed as 1.45 x 10 6 with 
U defined at station 1 and L = lm. The computations discussed below were based on 
the measured wall pressure distribution displayed in figure 6. Our curve fit by 
means of cubic splines approximated the measurements well, but we encountered prob- 
lems in obtaining smooth distributions of the pressure gradients, especially for the 
downstream region of the flow, where too few measurements have been carried out. 
Therefore using the plane-of-symmetry solution as starting conditions for the cross- 
wise integration (applied for case A4 of reference 4, see fig. 1) was not suffi- 
ciently accurate, but generated_large oscillations of all flow quantities. The 
profile of the normal velocity V at the first location off the plane of symmetry 
z = Az indicated a strong flow acceleration resulting in negative velocities V. 

This indicated that the profiles of 3W/3z calculated at z = 0 with 3 2 p/8z 2 and 
the crosswise pressure gradient prescribed at z = Az/2 were not compatible. Maybe 
smaller step sizes Az close to the plane of symmetry would have removed this problem. 
Instead we solved the boundary- layer equations with Hall's scheme in the vicinity of 
the plane of symmetry. Since in that region the near-wall turning angle of the flow 
was small, we shrank the computational domain with tan -1 (Az/Ax) = 3°. In. test runs 
we reduced the computational domain with tan -1 (Az/Ax) = 30° downstream of 
x = 350 mm and found negligible effects. The downstream side boundary of the com- 
putational domain was prescribed with z = 128 mm. 

The results obtained from calculations with the baseline turbulence model are 
indicated by solid lines ill figures 7 to 10. Up to station 4 the inner-layer pro- 
file of the mean velocity U-, figure 7, agreed well with the experimental data, while 
the outer-layer profile was underpredicted by about 10% at station 4. Further down- 
stream the near-wall flow was predicted to be larger than the measurements suggest- 
ing that for this part of the flow the turbulent viscosities were overestimated by 
equation (2). Indeed, at stations 5 and 6 the measured mixing length distributions 
could not be represented by £ = 0.41 y anymore but were much less. The mean veloc- 
ity component W was obtained close to the experimental data except for the outer 
layer of the profiles farthest downstream where the calculations were a few percent 
too large. The momentum thicknesses defined with the strearawise velocity profile in 
the direction of the outer-edge flow were found to be up to 10% higher than in the 
measurements. The calculated shape factor H was 5% larger than the experimental 
value at station 3 (H = 1.48), was close to the measurements at station 4 and 5 
(H = 1.47; H = 1.52), but about 7% smaller at station 6 (H = 1.56); At x = 500 mm 
and z = 50 mm a shape factor larger than 1.8 indicated incipient separation, which, 
in fact, occurred immediately downstream of this cross-section ( x = 5 08 mm and 
z » 44 mm). In figure 8 the calculated Reynolds shear stresses u m v m a nd v m w m are 
compared with the experimental data. In the outer layer the predicted UmVm correla- 
tion generally was too large. While at station 3 the peak value was only slightly 
higher than in the experiment, the maximum was too high by 20% at station 5 and 25% 
at station 6. The measured Reynolds stresses v m w m were not well represented by the 
computations; as for the van den Berg et al. experiment they were considerably over- 
estimated. However, the magnitude of this shear stress was much smaller than that 
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of the streamwise component (less than 1/4 up to station 4), and its effect on the 
mean flow development was expected to be fairly small. Additionally, the near-wall 
measurements tend to be smaller than zero. Such results are not compatible with 
mixing length theory and indicate severe experimental uncertainties. 

At stations 3 and 4 (x = 350 mm and 425 mm) the computed skin friction coeffi- 
cients as displayed in figure 9, were smaller than the Preston tube measurements by 
4% or 7%, respectively, while at station 5 (x = 475 mm) they were 10 percent larger. 
Since at station 6 the measured mean velocity profile defined in the direction of 
the near-wall flow, could not be described by the law of the wall anymore, the meas- 
urement is probably in error. In figure 9 we have also included the results of 
computations for the plane of symmetry, where two additional measurements were 
available. From the oil flow pattern of reference 2 we extracted several spanwise 
distributions of the wall turning angle, figure 10. We found good overall agreement 
even in the strong crossflow region except very close to separation (x = 500 mm and 
z < 65 mm) . 

In additional computations we investigated the sensitivity of the predictions to 
the prescribed turbulent length scale distribution. The experimental results did not 
indicate the outer layer similarity presumed by equation (2), but yielded £ e /6 = 0.07 
up to station 4 and Z /S « 0.06 further downstream. Approximating the measurements 
by a fixed length scale Z = 2.4 mm we obtained the results displayed_by broken 
lines in figures 7 to 10. The profiles of the mean velocity component U indicated 
considerably smaller velocities in the near-wall region. Correspondingly, the flow 
separ ated ahead of station 6 already. The magnitude of the Reynolds shear stress 
u m v m was considerably reduced, and the profiles were close to the experimental data. 
The crosswise shear stress only indicated a minor sensitivity. Corresponding to the 
mean-velocity profiles, the wall shear stresses were predicted lower than in the 
baseline computation (approximately by 0.00015 for x > 350 mm). The wall turning 
angle was close to the previous run up to x = 425 mm. At x = 475 mm the results 
were 5 larger at z max , but were 25° larger at z m ^ n and indicated incipient separa- 
tion. 


As investigated by further test calculations (not shown) , we found that the 
predicted near-wall flow quantities in the region close to separation significantly 
depended on the prescribed pressure distribution, which was measured to be larger 
at the wall than at the outer edge for x > 350 mm. Using the latter data together 
with the baseline turbulence model yielded a constant spanwise distribution of B w at 
x = 500 mm, and the flow separated further downstream compared to the results 
obtained with the wall pressure distribution. 

Generally, our numerical simulation of Dechow's experiment yielded reasonable 
overall agreement with the measured mean-flow profiles and wall quantities even in 
the strong crossflow region with turning angles at the wall of about 30° . However, 
similar to the experiment of van den Berg et al., the present one also indicated 
upstream pressure disturbances caused by flow separation on the test surface and 
probably enhanced by displacement effects of the flow at the upper-wall of the 
channel. Again, pressure gradients in the direction normal to the wall were present 
and thus impaired comparison calculations carried out by means of the boundary-layer 
equations. 
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MUller’s Experiment (Ref. 3) 


In this experiment a 3DTBL on a flat plate was generated by laterally deflecting 
the initial flow, figure 11. The static pressure distribution measured at the outer 
edge of the boundary layer is displayed in figure 12. We obtained a good representa- 
tion of the pressure field except for z > 500 mm by approximating the measurements 
by polynomial functions of fourth order. At five z-coordinates the distribution of 
Cp(x) was evaluated and then, at arbitrarily prescribed x-coordinates , the crosswise 
distributions Cp(z) were computed. To save computer time and because we did not 
encounter problems in calculating two-dimensional flows or those with moderate cross- 
flow, we did not start the numerical simulations at x = 0 but at x = 140 mm 
instead. At each crosswise grid point the smoothed velocity profile of station B5 
was normalized by the magnitude of the measured outer-edge flow and was assumed to 
be collateral to the outer-edge velocity vector. With this starting procedure we 
obtained velocity profiles at x = 200 mm, which differed less than 3% from the meas- 
urements, see, for example, profile B5 in figure 13. The lateral boundaries of the 
computational domain were defined by z max - 600 mm and by a wall streamline start- 
ing at x = 140 mm and z = 84 mm, while x max was prescribed with 650 mm. 

In figures 13 to 16 the measured downstream development of the flow is compared 
with the computed one obtained with the baseline turbulence model For all stations 
with z = 500 (row 5) the profiles of the mean velocity component U were in close 
agreement with the measurements, figure 13, while the peaks in the profiles of the 
W component were underpredicted by about 10%. Close to the separation line, at sta- 
tions D2 and E3, the agreement was noticeably worse. For z = 500 mm the momentum 
thicknesses evaluated for the streamwise profiles were found to be about 6% smaller 
than the experimental values. Since the displacement thicknesses were also under- 
predicted, the shape factors did not deviate more than 2% from the measurements 
(H = 1.46 at station E5 and H = 1.48 at station F5) . 


The computed and measured Reynolds stresses u m v m and v m wm are compared in figure 
14. Generally, good overall agreement was obtained for the streamwise component. 

The deviations for y < 5 mm at the upstream stations were prob ably caused by experi- 
mental inaccuracies. The profiles of the measured shear stress v m w m developed peaks 
at y ** 7 mm which were not obtained from the computations and were underestimated by 
as much as 30% at station E5 and 50% at station E3. However, since in the near-wall 
region the gradients of the computed crossflow velocities were underpredicted also, 
we could not expect closer agreement of measurements and predictions. In the outer 
layer the comparison yielded satisfactory results. 

The computed skin friction coefficients based on the total wall shear stresses 
are displayed in figure 15. For x > 200 mm the agreement was good except for 
z > 500 mm, where the outer-edge pressure distribution had been prescribed too high. 
At x = 600 mm the differences between measurements and computations increased with 
decreasing distance from separation and were as high as 15% at z m in* In figure 16 
the directions of the resultant wall shear stresses, computed at the measuring sta- 
tions, are compared with the wall flow visualization. Excluding the results at 
stations close to separation and those for z > 500 mm, the agreement was good, 
particularly since the crossflow angles 3w were as high as 40 . Generally, the 
results obtained from the numerical simulation of this experiment indicated that 
except for the region close to separation the downstream development of this flow 
was predicted well. 
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CONCLUDING REMARKS 


Numerical simulations of the 3DTBL-experiments of van den Berg and Elsenaar 
(ref. 1), Dechow (ref. 2) and MUller (ref. 3) were carried out in order to test the 
applicability of an algebraic mixing length turbulence model in flows with strong 
crossflow. Additional to the baseline integrations of the boundary- layer equations 
we also investigated the sensitivity of the predictions to modification of the tur- 
bulence model or to the prescribed outer-edge boundary conditions. Generally, for 
the regions with moderate skewing (g w < 15°) the flow-field calculations agreed well 
with the experimental data. For the highly three-dimensional parts of the flows, 
however, we found different degrees of agreement for each case. 

The computed downstream development of the flow measured by van den Berg and 
Elsenaar failed to simulate the measurements. Separation was only obtained if the 
pressure distribution prescribed was larger than the one measured at the wall and 
simultaneously the magnitude of the eddy viscosities was reduced compared to the 
baseline model; this was also true in the computations of Elsenaar et al. (ref. 10). 
Under these conditions the results proved to be extremely sensitive to either modifi- 
cation; for example, changing the pressure distribution by a few percent yielded 
large changes in the computed flow field. We feel that the problems encountered were 
these conditions the results proved to be extremely sensitive to either modification; 
for example, changing the pressure distribution by a few percent yielded large 
changes in the computed flow field. We feel that the problems encountered were 
mainly a result of the breakdown of the boundary-layer concept in the downstream part 
of the flow. Though in the experiment the diffuser cross-section approximately 
varied linearly throughout the flow field, the longitudinal pressure gradients 
rapidly decreased beyond station 5. This upstream pressure effect was caused by 
separation and the flow blockage at the downstream end of the diffuser. Additionally, 
pressure variations normal to the wall probably influenced the development of this 
flow. For resolving these elliptic flow-field effects, predicting the location of 
separation and ultimately testing turbulence models, a Navier-Stokes prediction 
method would be necessary. 

Our numerical simulation of Dechow f s experiment yielded reasonable overall 
agreement with the measured mean-flow profiles and wall quantities even in the strong 
crossflow region with turning angles of about 30° between wall- and outer-edge flow 
direction. However, the flow field close to separation was found to be sensitive to 
the pressure distribution prescribed, as well as to the turbulence model used. In the 
three-dimensional flow, the measurements indicated normal pressure gradients which 
impaired the comparison with solutions of the boundary-layer equations. Additionally, 
in the downstream part of the flow similarity laws as developed for two-dimensional 
mean flow and turbulence fields were not valid anymore. Under these conditions we 
could not extract detailed information about turbulence modeling. 

In MUller ? s experiment, which was performed in a wind tunnel with an open test 
section of 1 m 0, most measuring stations were located at sufficiently large dis- 
tances away from separation in order to ensure negligible normal pressure gradients. 
The results of the numerical simulation for this flow indicated a minor dependency on 
the accuracy of the prescribed pressure distribution and enabled us to study the 
applicability of the turbulence model used. Except for the region closest to separa- 
tion, the flow development, computed with a baseline turbulence model, compared 
favorably with the measurements even at the measuring stations farthest downstream 
with large skewing (3 W > 40°). In contrast to the other calculations, the results 
were found to be insensitive to the outer-layer mixing length and suggested the 
applicability of an algebraic eddy-viscosity model as developed for two-dimensional 
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boundary layers. In order to fully understand the results obtained, further investi- 
gations with more sophisticated turbulence models will be necessary, which take into 
account the experimentally observed anisotropy of the turbulence field or the normal 
Reynolds stress diffusion, respectively. 

The present study again demonstrates the lack of well-documented experiments 
which are suitable to guide both the development of computational methods for 3DTBLs 
and closure assumptions as well. Because pressure forces play a dominant role in 
flows like those investigated here, future experiments should accurately map out the 
pressure distribution and also ensure that the boundary layer concept is valid. 

Based on such experiments we expect progress in turbulence modeling for 3DTBLs as 
well as in the development of prediction methods for engineering-type applications. 
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Figure 4. - Comparison of measured and computed mean velocities. 








Figure 7. - Comparison of measured and computed mean velocities. 
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Figure 9. 


- Comparison of measured and computed skin friction coefficients 
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Figure 10. - Comparison of measured and 

computed wall turning angles. 
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Figure 11. - Experimental setup of reference 3. 
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Figure 12. - Measured and smoothed pressure distribution. 
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Figure 14. - Comparison of measured and computed Reynolds shear stresses. 
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Figure 16. Comparison of oil flow pattern with wall— streamline directions computed at measuring stations 
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